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Abstract 

We present a generalization of the Sutherland's multicomponent model. 

Our extension includes both the ferromagnetic and the antiferromagnetic t-J 

model for any value of the exchange coupling J and the hopping parameter t. 

We prove rigorously that for one dimensional chains the ground-state of the 

generalized model is non-degenerate. As a consequence, the ordering of energy 

levels of the antiferromagnetic t-J model is determined. Our result rigorously 

proves and extends the analysis carried out by Sutherland in establishing the 

phase diagram of the model as a function of the number of components. 
74.20.-z, 71.20.Ad, 75.10.Jm 
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The t-J Hamiltonian is certainly one of the most attractive model to explain high tem- 
perature superconductivity (HTc) without any phonon pairing mechanism Despite its 
simplicity and the huge amount of literature presented starting soon after the discovery 
of HTc, the determination of the full phase diagram, especially in two dimensions, is still 
controversial. 

Recently, after a work of Bares et al. it was immediately realized that the t-J model 
at the supersymmetric points J = ±2t is a particular case of a multicomponent model 
introduced previously by Sutherland || in his investigation of models solvable by Bethe 
Ansatz. In the present work we generalize this mapping as to include the general and 
physically more interesting J ^ ±2 1 case. Our approach allows the determination of some 
rigorous and novel result about the t-J Hamiltonian and other related models. We also 
discuss in some detail the relation between the multicomponent model and the usual form 
of the t-J Hamiltonian, trying to clarify some delicate point concerning the interpretation 
of Sutherland's work and his notations. 

Consider a hypercubic lattice of N sites and place N objects, or particles, on it exactly 
one to a site (hard-core condition). The particles are of z different species, of which x are 
assumed bosonic and y fermionic. We define the multicomponent t-J model 

tf = -*£^-££J#> (!) 

<ij> Z <ij> 

where the sum is extended to nearest-neighbor (NN) sites, and where Pfj and permute 
particles occupying sites i,j if these have opposite or same statistics, respectively, otherwise 
they give zero. The operators in Eq. ([!]) are thus defined via the decomposition P^ = 
+ P^ of the operator P it j permuting whatever pair of particles occupies sites so 
that for J = It Eq. (|1|) reduces - apart for an overall constant - to the Sutherland's 
multicomponent Hamiltonian H$ = — J2<ij> P,j- We notice that H is equivalent to H$ also 
for the case of a purely bosonic (y = 0) or fermionic (x = 0) system, as for in this case 
P/j = 0. The Hamiltonian ([!]) is better interpreted as a model for a single boson B and 
a single fermion F having x and y different colour states, respectively. According to this 
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simple interpretation, and following the notation of Sutherland, we refer to Hs as the B x F y 
system and similarly we refer to Eq. (|IJ) as the B x Fj system. In the latter case we employ 
the subscript J to stress that the Hamiltonian ([]]) carries a free parameter (i.e., the ratio 
J/t) not present in the supersymmetric model Hs- Henceforth we consider only bipartite 
lattices. In this case the sign of the hopping parameter in Eq. (1) is irrelevant, because 
on a bipartite lattice it can always be changed by means of a unitary transformation [[|. 
Henceforth, we assume t > 0. More interestingly, also the sign of the coupling J can be 
fixed without loss of generality, because the B x Fj and B y F x j systems are equivalent ||. 
Henceforth in Eq. (|l|) we set J > 0. Using this sign convention the BF 2 and B 2 Fj systems 
are (apart from a shift in the zero of the energy) the usual t-J Hamiltonian for J > and 
J < 0, respectively, whereas the F 2 and B 2 systems correspond to the antiferromagnetic and 
ferromagnetic Heisenberg models. The BF system is instead a simple free spinless fermion 
Hamiltonian. One can also easily check that the B 2 F 2 system is the EKS model introduced 
by Efiler et al. || in their investigation of ^-pairing superconductivity, and that the B 2 F 2 =0 
system corresponds to a bond-charge Hubbard model as discussed by Arrachea and Aligia 

To avoid a too cumbersome notation, we employ a single colour index a = 1, 2, .., y, y + 
1, ..,z, with first y and last x = z — y labels referring to fermion and boson colours, re- 
spectively. To distinguish the statistics we then set £ a = 1 (—1) for bosons (fermions). 
Henceforth we denote N F = YZ>=i ^ a an d N B = J2a= y +i N a the total number of fermions 
and of bosons, respectively, where N a is the number of particles of colour a. We remark that 
each number N a is individually conserved because in Eq. ([!]) there only enter permutators. 
An immediate and important consequence which will be used later is that: 

(A) Each B n F" l J l system with n < x and/or m < y is a particular subspace of the B x Fj 
system Hilbert space, 

which is obtained by setting to zero the occupation numbers of (x — n) boson and (y — m) 
fermion colours. 
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To describe the model ([!]) in second quantization we employ a simple oscillator rep- 
resentation where y ai (f a ,i) creates (destroys) a particle of colour a on site i out of a 
unphysical reference vacuum |0), and accordingly f^, f a ^ are assumed as fermion (boson) 
operators for £ a = — 1 (£ Q = 1). The hard-core condition leads to the set of local constraints 
Sa=i flifa,i = 1> Vi, so that a vector of the physical Fock space is given by a product of N 
distinct creation operators acting on the vacuum. On any of such Fock vectors the operator 
Pij acts by interchanging the site indices 

i', r ..f;,. r .,i\ l ...o) ..../;;.,..,/•,,...()). (2) 

(i.e., it acts like a "Majorana" permutator). In the physical Fock space we thus have 

= fa,jf/3,if/3,jfa,i 

= fa,jf/3,if/3,jfa,i (3) 

We notice that if one forces to represent Pij by means of the Dirac permutator P/? in- 
terchanging the colour labels (i.e., Pijfa,ifp,j\ty = flifa,j\ty)i then an extra minus sign is 
involved for the exchanges of two fermions (i.e., P^j = —Pfj when permuting two fermions, 
and Pi j = P® otherwise). Due to this property, Pij is also defined in literature as a graded 
permutator PJ^]. Henceforth we always refer to Eq. (0), because this is the definition con- 
sistent with Sutherland's notations. 

In the oscillator representation Hubbard operators have the simple bilinear expression 
= faifp,i- As is well known [jj, the total operators T a @ = J^Li f° r both fermionic 
(£<* = 6/3 = — 1) an d bosonic (£ a = ^ = 1) labels are the generators of two commuting 
SU(y) and SU(x) algebras, respectively. The diagonal operators of the algebras are the 
numbers N a , i.e., T a a = N a , whereas the off-diagonal operators T a p, for a ^ (3, represent 
Cartan raising-lowering step operators, as for they change by ±1 the numbers N a ,N^. As 
H commutes with both the algebras, 

[H,T a p} = for £a = £/3, (4) 
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it follows that for each fixed value of Np (hence, of Nb = N — Np) the spectrum of H can 
be classified in multiplets of SU(?/)(8>SU(x) ||. Our task is to characterize the sector of the 
ground-state and, more in general, the ordering of the energy levels of the model (|I|) in one 
dimension. To this aim, we impose open boundary conditions, because this choice greatly 
simplifies the discussion. 

• We consider the sector {N a } with definite values of the numbers N a . The main result 
of this paper is the proof that: 

(B) For non-vanishing t and J, in any of the possible sectors {N a } the ground-state of the 
Hamiltonian (fjj on an open chain is non-degenerate. 

Denoting ii,..,i N i the positions of the particles of colour a — 1, then i N i + i, ..,i N i +N 2 
those of the particles of colour a = 2, and so on, a basis {|a)} of the sector is given the 
N\/(J\a=i N a \) distinct sequences 

that we order according to the fundamental regions 



i\ < %2 < ... < ijyi 
«jv!+i < *at 1 +2 < ••• < im+N 2 



(6) 



iN-N z +l < ^N-l < ••• < 

where all operators of same colour are grouped together, and where colour labels and, 
separately, coordinates of particles of same colour are ordered in ascending order. Because 
of the open boundary conditions, one can easily see that in the chosen basis the off-diagonal 
matrix elements H a ^ = (a\H\b), a ^ b, of the Hamiltonian (0) are non-positive. From this 
property, by evaluating the expectation value of the energy over a normalized state-vector 
1^) = T,ala\a) one obtains || 

^ = E^a, 6 7:76>E^lTa||T6|, (7) 

a, b a, b 



so that a ground-state can be chosen with non-negative coefficients, 7 a = |7 a |. 

The non-degeneracy of the ground-state is then proved by showing that for the ground- 
state the coefficient of each element of the basis is non-zero, 7 a = |7 a | 7^ 0, Va, as for 
no other ground-state with the same restrictions on the 7 a 's can satisfy the orthogonality 
requirement. According to Perron- Frobenius theorem, the above more restrictive result 
holds provided two arbitrary states \a), \b) are linked by the Hamiltonian ||10|| , namely one 



can find an integer Q such that (b\H®\a) 7^ 0. The Hamiltonian H satisfies this property 
whenever t and J are both different from zero. In fact, the site indices {jk} of the state 
1^) = fi,jr--fz,j N \fy define a permutation P of those, {ik}, of the state \a). As is well 
known, one can always resolve the total permutation P as a product of NN permutations 
P = rifcLi Pi h ,i k +i, i- e -? \b) = P\a). The Hamiltonian (P is the sum H = — YnLiP^i+i of 
NN operators V iji+1 = tP* i+l + (J/2)P$ +1 . By definition, the expansion of H® contains 
the product V = E[fc=i wnere {h} is exactly the sequence of site indices entering the 

decomposition of P. If t and J are both non-zero one has V\a) = uP\a), where v - if Q is 
assumed as the minimum integer allowing the decomposition - is a strictly positive constant. 
Hence, we obtain |(&|if^|a)| > u, because all other sequences of permutations contained in 
the expansion of H® that in case link the states \a) and |6) contribute to the expectation 
value of (b\H®\a) with the same sign. This concludes the proof of non- degeneracy of the 
ground-state for t 7^ and J 7^ 0. 

• The uniqueness of the ground-state (henceforth we avoid the case t = or J = 0) is 
valid in any sector with fixed colour numbers, but of course one can have many degeneracies 
among distinct sectors. Indeed, the symmetries (P allow a simple proof that: 

(C) For an open chain, the ground-state energies of the B x Fj and of the BFj system are 
degenerate. 

In fact, from property (A) the ground-state |^o) of the B x Fj system in the N y+l = N y+2 = 
... = N z _i = 0, N z = Nb sector is the ground-state of the BFj system. This state 
is a maximal weight for the algebra SU(x) annihilated by all bosonic raising operators, 
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F z a\i>o) = 0, for a < z, £ a = 1, so that by applying on it bosonic lowering operators Y a z one 
can obtain a degenerate eigenstate of any sector of the B x Fj system with arbitrary bosonic 
occupation numbers N a but fixed Nb- These kind of states satisfy the condition 7 a = |7 a |, 
and from the uniqueness property we conclude that each of them must be the ground-state of 
its sector. This concludes the proof. By taking the thermodynamic limit one can then easily 



show that the energy per site does not depend on the boundary conditions |TT[, confirming 
and extending to the general J ^ 2t case the Sutherland conclusion || that the B x F y and 
the BF y systems have not only the same ground-state energy per site but also the same 
ground-state (both equivalent by the SU(a;) symmetry). This statement however holds only 
in one dimension, and cannot be extended to higher dimensions, where counter-examples 
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• The uniqueness of the ground-state also implies that all its quantum numbers are 
well defined. As usual in this case, these numbers can be identified by selecting some simple 
reference state = J2 a Pa\ a ) (e.g., ground state of a simple reference Hamiltonian) with 
non negative amplitudes p a = \p a \ and for which the quantum numbers are known. By the 
non vanishing value of the scalar product (</>|^) = J2a \Pala\ h follows that the ground-state 
must belong to the same sector. 

The discussion of the ordering of energy levels and their associated quantum numbers 
in its generality is quite involved. For the sake of clarity, we now specialize to the simpler 
but physically most interesting case of the B 2 F 2 system. For convenience we rename the 
Hubbard operators 

t+,i — r i j2) T-,i = r^i, T x>i = — (r^ ^ - v 2 2 ), (8) 
= r? 4 , l_j = r^ 3 , L Zti = -(r? 3 — r^ 4 ), (9) 



and denote Tj = (T Xji , T y>i , T Z)i ), where T±^ = T Xji ±iT yji , and similarly for Lj. We also denote 
T = Y^iLi Ti and L = YliLi Li, which thus are the fermion and boson total spin operators, 
respectively. The eigenstates are classified according to the quantum numbers Np, T z , L z , 
T, L, where T(T + 1) and L(L + 1) are the eigenvalues of T 2 = T ■ T and L 2 = L ■ L, 



respectively. 

Result (C) and Eq. ([]) imply that the ground-state of the B 2 F] system is characterized 
by maximum value of the boson total spin, L = \L Z \ = Nb/2. To discuss the value of the 
fermion total spin in the ground-state, we follow the approach of Lieb and Mattis || for the 
case of the Heisenberg model. After division of the chain into two sublattices A and B, we 
define the reference Hamiltonian 

H = 2T A -f B , (10) 



where T p = ^ i£p Ti (p = A,B). The Hamiltonian (|10|) commutes with local fermion n{ = 
Y\ 1 + T 2 2 and boson n\ = Ff 3 + Ff A occupation numbers. Denoting nf p = J2ie P n { the 
number of fermions on sublattice p, the lowest eigenvalue E(T, L) of H for each given T and 
arbitrary L is easily found 

E(T, L) = T(T + 1) - - £ nf(nl + 2), (11) 
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^ p=A,B 



where the allowed values of T are 



n{-nt\ m N F 

< T < — . (12) 



2 ~ ~ 2 

In the basis (^) the sign of the off-diagonal matrix elements of Eq. fllCf ) is not definite. 
However, to select a good reference state it is sufficient to consider the subset Q of the 
basis vectors where the fermions occupy the first iV^ sites of the chain. In fact, it is not 
difficult to check that in the invariant subspace Q the off-diagonal matrix elements of H 
are non-positive, because in Eq. fllO"D permutations take place only between sites of different 
sublattices. The previous analysis (0) implies therefore that the amplitudes p a of the lowest 
energy eigenstate can be chosen non-negative. We conclude that the ground-state of H has 
the same quantum numbers of the reference state. In the subspace Q one has 

f N F f N F 
n J A = ——, n B = —— tor JMp even, 



and by minimizing Eq. ( pTTj) one finds that the ground-state of the B 2 F] system has T = 
(for Np even) or T = 1/2 (for Np odd). More in general, following Ref. ||, from Eq. ( [TT]) 
we can order the energy levels E(T, L) of the B 2 F] system in function of the fermion total 
spin quantum number 

E(T, L) < E(T + 1, L). (14) 

• Eq. (|i~4f) and the result (A) allow the determination of the ordering of the energy levels 
Ej for the antiferromagnetic t-J Hamiltonian. 

The correspondence between the BF] system and the t-J model for J > is explicitly 
derived upon identification of spin-up and spin-down electrons with the two fermion colours 
and of the hole with the boson colour a = 3. Accordingly, we have T = S and Nb = 
Nh, where S and Nh are the physical spin and hole number, respectively. Because BF] 
is the subspace of B 2 F] with L = \L Z \ = Nb/2, from Eq. ([14]) we derive the ordering 
E J>0 (S, Nh) < E J>0 (S + l,N h ). In particular, we find that for J > the ground-state of 
the t-J model on an open chain has minimum total spin S and is non degenerate (modulo 
the trivial twofold S z = ±1/2 spin degeneracy for Nh odd). 

In a similar fashion, the correspondence between the B 2 Fj and the t-J model for J < 
is obtained by identifying spin-up and spin-down electrons with the two boson colours and 
the hole with the fermion colour a = 1. In this case we have L = S and Np = Nh, so 
that by applying result (C) we find that for J < the ground-state of the t-J model on an 
open chain is [modulo the trivial (N — Nh + l)-fold spin degeneracy] a nondegenerate fully 
polarized ferromagnet. 

Concerning the general case, we emphasize here that the degeneracy among sectors with 
different number of bosonic colours but with the same Nb is due to the fact that the quadratic 
operator 

a,/3>y 

of the bosonic algebra assumes its maximum value in the ground-state. Concerning the 
fermionic SU(y) sector for generic y > 2, symmetry considerations, the case y = 2, as 
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well as the exact solution at the supersymmetric point ||, suggest that the ground-state of 
the Hamiltonian (|1|) should have minimum value of the quadratic operator Cf [namely, the 
operator as defined in Eq. ( |i5|) but for fermions, i.e., a, (3 < y}. In particular, for values of the 
fermion number Np = yn, with n an integer, we conjecture that the absolute ground-state 
of the multicomponent model ([]]) necessarily belongs to the sector with N l = ... = N y = n 
for any value of the coupling J > 0. 

We have presented a multicomponent generalization of the t-J model, extending in a 
natural way the well known Sutherland's Hamiltonian corresponding to the supersymmetric 
J = ±2t cases. In one dimension we have determined that doping does not alter the 
ground-state spin sector of the Heisenberg Hamiltonian, nor the ordering of energy levels of 
the antiferromagnetic model. These statements have been derived by using Perron-Frobenius 
theorem, i.e., by exploiting the irrelevance between fermion and boson statistics, typical for 
one dimensional systems with NN interactions. Indeed, the off-diagonal matrix elements 
of the Hamiltonian ([]]) can be made non-positive by a proper choice of the fermion-boson 
basis. This property is clearly lacking in higher dimensions, and consequently the equivalence 
between the B x Fj and the BFj systems, valid in one dimension, is not rigorous in higher 
dimensions - unless for the trivial case without fermions. In fact, contrary to what previously 
stated HI, it is easy to find a counter-example, where the ground-state of the B 2 Fj system 
(i.e., the ferromagnetic t-J model) does not belong to the sector of maximum spin BFj, 



[i.e., to the Nagaoka sector) but to the singlet sector |12 
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